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ABSTRACT 

Signals  generated  by  moving  traeked  vehieles  ean  be  signifieantly  affeeted  by  near-surfaee 
materials  ineluding  pavements,  eompaeted  soils,  shallow  bed  roek,  and  the  water  table.  The 
ways  that  seismie  signals  depend  on  these  surfaee  layers  are  examined  by  analyzing 
impulsive  point  forees  applied  at  the  earth-air  interfaee.  Fields  are  ealeulated  using  a  plane 
layered  frequeney  domain  wavenumber  model.  Synthetie  time  domain  seismograms  are 
generated  by  eonvolving  the  layered  earth  impulse  responses  over  a  band  of  frequeneies  with 
an  appropriate  souree  foree  funetion.  The  effeet  of  eonstant  Q  attenuation  is  also  ineluded. 

Our  results  are  related  to  seismie  signals  from  traeked  vehieles  by  noting  that  the  spatially 
distributed  nature  of  the  traek-ground  foreing  at  any  instant  in  time,  ean  be  obtained  by 
superposition  of  multiple  point  souree  results.  In  generating  earth  impulse  responses,  we  use 
surfaee  layering  that  is  eommonly  eneountered  moving  vehiele  problems.  In  eaeh  layer  ease 
eonsidered,  we  systematieally  vary  layer  thiekness  to  demonstrate  how  thiek  a  partieular  layer 
must  be  to  generate  a  signifieant  alteration  in  the  eharaeter  of  the  seismie  signals.  It  is  found 
that  surfaee  layers  less  than  about  a  meter  thiek  do  not  have  a  major  effeet  on  observed 
waveforms  in  the  frequeney  bands  of  interest.  However,  for  near  surfaee  layers  greater  than 
1.5  m,  eomparisons  between  high  veloeity  and  low  veloeity  materials  ean  have  amplitude 
differenees  as  great  as  20  dB  at  ranges  of  100  m.  Waveforms  from  these  two  models  also 
have  substantively  different  surfaee  waves  modal  eharaeteristies. 

1.  INTRODUCTION 

The  seismie  signals  reeorded  from  moving  vehieles  ean  be  used  to  traek,  and  identify  them.  However, 
for  seismie  signals  to  aehieve  their  full  potential  in  fieldable  battlefield  systems,  the  wide  variety  of 
geologie  effeets  on  signal  properties  must  be  eonsidered.  The  upper  10  m  of  the  geologie  materials  ean 
have  a  partieularly  pronouneed  effeet  on  propagating  seismie  surfaee  waves.  The  effeets  inelude  large 
variations  in  signal  amplitude,  and  waveform  dispersion.  The  eharaeter  of  the  reeorded  surfaee  wave 
signals  ean  be  thought  of  as  a  eonvolution  of  the  foree  the  vehiele  applies  to  the  earth  and  the  response  of 
the  earth  to  the  foreing.  In  this  paper  we  systematieally  apply  a  layered  media  seismie  propagation  model 
to  several  eommon  geologie  seenarios  and  analyze  the  effeets  of  layer  properties  on  the  seismie 
waveforms. 
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2.  METHOD 

Seismograms  are  eomputed  in  two  steps.  First  at  a  number  of  frequeneies,  the  Green’s  Funetion  is 
ealeulated  for  the  response  of  a  point  souree  in  an  infinite  media  by  using  a  wavenumber  integration 
method.  The  geometry  of  this  ealeulation  is  shown  in  Figure  1.  The  foree  time  funetion  of  the  souree  is 
then  transformed  into  the  frequeney  domain.  Next,  the  Green’s  Funetion  and  the  transform  of  the  time 
foree  funetion  are  multiplied  together  in  the  frequeney  domain  (this  gives  a  eonvolution  in  the  time 
domain).  Finally,  the  result  is  transformed  baek  to  the  time  domain  to  generate  seismogram  waveforms. 

The  most  diffieult  part  of  the  eomputation  is  generating  the  frequeney  domain  Green’s  Funetion.  This  is 
done  using  the  program  GRPOINT  (Hisada,  1994;  Hisada,  1995a;  Hisada,  1995b).  This  program  gives 
the  response  for  direeted  point  forees  at  arbitrary  depth  in  a  layered  media;  the  depth  of  the  reeeiver  is 
also  arbitrary.  Both  vertieal  and  horizontal  direeted  point  forees  are  allowed.  All  three  eomponents  of  the 
ground  motion  are  ealeulated.  A  program  was  developed  to  take  the  output  of  GRPOINT  put  in  a  number 
of  different  foree  time  funetions,  and  eompute  seismograms. 

In  this  paper,  unless  indieated,  all  results  are  for  a  vertieal  point  foree  souree  (TX)  at  the  surfaee  with  a 
reeeiver  (RX)  on  the  surfaee.  The  seismogram  give  the  vertieal  veloeity  of  ground  partieles.  This  is 
analogous  to  the  signals  reeorded  by  vertieal  geophone.  The  results  are  denoted  as  wavenumber 
integration  (WI)  results.  The  foree  time  funetion  of  the  vertieal  point  foree  is  shown  in  Figure  2. 

3.  RESULTS 

To  verify  the  WI  ealeulation  a  eomparison  was  made  with  time  domain  finite  differenee  (TDFD)  method 
results  obtained  with  the  PTOP  program  (Keteham  et  al,  2000;  Hestholm  and  Ruud,  1998  ).  The  results  , 
shown  in  Figure  3,  show  the  two  methods  give  very  similar  waveforms,  and  thus  both  are  giving  aeeurate 
results. 

The  seismograms  for  the  model  L2,  shown  in  Figure  4,  are  typieal  for  a  low  veloeity  over  a  halfspaee. 

The  major  arrival  is  a  dispersed  fundamental  Rayleigh  wave  arrival.  The  seismograms  for  the  model  LI, 
shown  in  Figure  5.  In  this  model  the  near  surfaee  layer  has  a  very  large  s-wave  veloeity  eontrast  with  the 
half  spaee,  whieh  introduees  a  sealloped  appearanee  in  the  waveforms.  This  is  attributed  to  the 
interferenee  of  the  fundamental  Rayleigh  wave  mode  with  a  higher  order  Rayleigh  wave  mode.  The 
result  sharply  eontrasts  with  the  mueh  simpler  waveforms  seen  in  the  L2  model. 

The  next  model,  HI,  has  a  high  veloeity  layer  at  the  surfaee;  waveforms  are  shown  in  Figure  6.  Sueh  a 
model  eould  represent  a  frozen  layer  above  uneonsolidated  soil.  The  waveform  for  HI  show  two  distinet 
frequeney  bands.  We  hypothesize  that  the  early,  high  frequeney  arrival  is  energy  that  has  traveled  in  the 
layer,  while  the  later  low  frequeney  wave  is  the  fundamental  Rayleigh  wave  in  the  whole  layered  system. 

Figure  7  gives  a  eomparison  of  the  amplitudes  versus  distanee  for  the  three  models.  The  L2  model  has 
the  highest  amplitudes,  the  LI  model  is  next;  the  HI  model  gives  the  smallest  amplitudes.  At  a  range  of 
1 10  m  the  L2  model  amplitude  is  about  20  dB  above  the  HI. 

Using  the  low  veloeity,  LI,  geologie  model  we  explore  the  effeet  of  layer  thiekness.  The  results  are 
shown  in  Figure  8.  For  H  =  0  there  is  no  low  veloeity  layer  and  a  simple  Rayleigh  wave  pulse  on  the 
surfaee  of  a  high  veloeity  material  halfspaee  is  observed  at  about  .1  s.  The  effeet  of  the  layer  does  not 
alter  the  waveform  signifieantly  until  H  approaehes  1  m.  Beyond  this  thiekness,  a  slower  arrival  starts  to 
appear.  This  is  energy  that  is  traveling  in  the  layer  as  a  leaky  guided  mode.  The  sealloped  waveforms  are 
observed  from  interfering  Rayleigh  modes  are  observed  for  1  m  <  H  <  8  m.  At  H  =  8  m,  the  high  speed 


early  arrival  is  very  weak  eompared  to  the  Rayleigh  wave  arrival.  Finally  for  an  infinitely  thiek  layer  (  H 
=  oo  )  a  simple  Rayleigh  wave  arrival  is  seen. 

Using  the  HI  model,  the  effeet  of  high  veloeity  layer  thiekness  on  waveforms  is  shown  in  Figure  9.  For  H 
=  0  there  is  no  high  veloeity  layer  and  a  simple  Rayleigh  wave  pulse  on  the  surfaee  of  a  low  veloeity 
material  halfspaee  is  observed  at  about  .34  s.  The  effeet  of  the  layer  does  not  alter  the  waveform 
signifieantly  until  H  approaehes  0.75  m.  At  thiekness  larger  then  this,  a  higher  veloeity  arrival  begins  to 
appear.  This  is  energy  that  is  traveling  in  the  layer.  The  waveform  for  H  =  1  m,  shows  the  early  higher 
veloeity,  high  frequeney,  arrival  and  the  later  lower  veloeity  arrival.  Finally  for  an  infinitely  thiek  layer  ( 
H  =  oo  )  a  simple  Rayleigh  wave  arrival,  for  the  high  veloeity  material,  is  seen. 

4.  CONCLUSIONS 

We  have  eompared  our  wavenumber  integration  results  with  ERDC-CRREL’s  finite-differenee  models 
and  demonstrated  exeellent  amplitude,  arrival  time  and  frequeney  eontent  agreement.  This  demonstrates 
that  both  methods  are  eable  of  generating  eorreet  results.  In  general,  our  systematie  analysis  of  simple 
layered  models  shows  that  strong  veloeity  variations  of  near  surfaee  layers  will  not  effeet  signal 
eharaeteristies  until  the  layer  thiekness  exeeeds  1  m.  This  is  partieularly  true  of  high  veloeity  materials. 
Low  veloeity  geologie  materials  overlaying  high  veloeity  materials  have  the  greatest  impaet  on  signal 
properties  by  supporting  multiple  high  amplitude  Raleigh  wave  modes.  In  this  ease,  waveforms  show 
extended  (in  time)  group  envelopes  with  eomplieated  internal  interferenee  eharaeteristies.  Lastly  we  note 
that  at  ranges  of  1 10  m  the  signal  amplitudes  ean  differ  by  as  mueh  as  20  dB  between  low  veloeity  layers 
and  high  veloeity  layers. 
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Figure  1.  Geometry  for  point  source  in  a  layered  medium.  Force  can  be  in  arbitrary 
direction,  r  is  the  range  to  the  RX. 
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Figure  2.  Force  time  function  and  it’s  power  spectra. 
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Figure  3.  Comparison  of  waveforms  from  PTOP  (FDTD)  and 
wavenumber  integration. 
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Figure  4.  Waveforms  at  varying  ranges  for  the  L2  model.  L2  model  has  a  small  velocity 
contrast  between  the  layer  and  the  halfspace.  H  =  1.5  m. 
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Figure  5  Waveforms  at  varying  ranges  for  the  LI  model.  LI  model  has  a  large  velocity 
contrast  between  the  layer  and  the  halfspace.  The  s-wave  contrast  is  more  than  a  factor  of 
four.  H  =  1.5  m. 
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Figure  6.  Waveforms  at  varying  ranges  for  the  HI  model.  H  =  1.5  m. 


Figure  7.  Amplitude  versus  distance  for  three  models.  H  =  1.5  m. 
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Figure  8.  Waveforms  for  varying  H,  for  the  LI  model.  H  given  on  figure. 
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Figure  9.  Waveforms  for  varying  H,  for  the  HI  model.  H  given  on  figure. 


